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Abstract 

Based on Friedberg and Lee's geometric picture by which the tribimaximal Pontecorvo-Maki- 
Nakawaga-Sakata leptonic mixing matrix is constructed, namely, corresponding mixing angles 
correspond to the geometric angles among the sides of a cube. We suggest that the three realistic 
mixing angles, which slightly deviate from the values determined for the cube, are due to a viable 
deformation from the perfectly cubic shape. Taking the best-fitted results of #12 and #23 as inputs, 
we determine the central value of sin 2 26\z should be 0.0238, with a relatively large error tolerance; 
this value lies in the range of measurement precision of the Daya Bay experiment and is consistent 
with recent results from the T2K Collaboration. 

PACS: 14.60.Pq Neutrino mass and mixing 
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I. INTRODUCTION 



Neutrino oscillation observations have revealed evidence that neutrinos are massive. Neu- 
trinos are produced via weak interaction as flavor eigenstates i/f = (z/ e , z/ M , v T ) and can be 
written in the mass basis v m = (ui, z/ 2 , z/ 3 ), which are really the physical states. These two 
bases are related by a unitary matrix U u , i.e., Vf = U v Vm. The mixing in the lepton sector is 
named as the Pontecorvo [jJ-Maki-Nakawaga-Sakata Q (PMNS) matrix, which can account 



for the currently available data on the observation of so 
and the reactor and accelerator neutrino experiments 
charged currents are 



ar, atmospheric neutrino oscillations 
"I 

3J. In the standard model, the weak 
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where i,j = 1,2,3 and correspond to physical particles. The mixing matrix 
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is a 3 x 3 unitary matrix and can be parameterized by three mixing angels #12, #23, and #13, 
and one CP phase S 
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where = cos^j, Sij = sin^-. If neutrinos are Majorana particles, there would be an ad- 
ditional diagonal matrix diag(e 4ai//2 , e ta2 ^ 2 , 1) multiplied to the above Upmns matrix, which 
is not relevant for neutrino oscillations. The parametrization Eq. ([3]) can be rewritten as a 
product of three rotations Rij in the ij plane through angles 6^ and a diagonal CP phase 
matrix Us=diag(e l5 ^ 2 , 1, e~ l8 ^ 2 ), 
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There have been numerous phenomenological Ansatze for the entries of Upmmsj for ex- 
ample, the democratic {4], the bimaximal [5|, and the tribimaximal Ansatze 6|. Among 
them, the tribimaximal mixing is closer to the experimentally observed mixing patterns, 
and the matrix is given by 



Utri 



tribi 



( i/Vs ^ 

-l/y/6 1/V3 l/y/2 

\y 1/V6 -1/V3 1/V2 j 



(6) 



which suggests 6>i 2 = sin -1 (l/\/3), #23 = 7r/4, #13 = 0. As noted, the CP phase e tS is 
always associated with S13 [Eq. (J3J)]; thus, null 6*13 would imply that one cannot observe 
CP violation at lepton sector in the framework of the standard model even though 5 is 
not zero. Obviously, there is no priori that the CP violation should appear at the lepton 
sector, but only nonzero #13 can intrigue an enthusiasm to explore CP violation at the lepton 
sector. Once the #13 is determined to be nonzero as the T2K experiment and our theoretical 
prediction made in this work suggest, the next step would be searching for CP violation at 
the lepton sector. 

Indeed, the tribimaximal symmetry is well manifested by the data. A rigorous symmetry 
would demand #13 to be zero; however, it is not the whole story because this elegant sym- 
metry is to be broken, and a nonzero #13 is expected. The question is if it is not zero, what 
is its size, which is the main concern of the recent studies. 

The unbroken tribimaximal matrix Eq. ([6]) can be further written as a sequential product 
of two independent rotations on 12 and 23 planes: 
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-l/y/3 2/y/E 

V V 

and -R13 becomes a 3 x 3 unit matrix, i.e., Utriu = -R23(7r/4)-Ri2(sin _1 (l/v / 3))- Friedberg 
and Lee 7| propose a geometrical interpretation for the tribimaximal symmetry as shown 
in Fig. [TJ For readers' convenience, let us briefly introduce Friedberg and Lee's geometrical 
model and their conventions The charged leptons in the basis L c = (S e , S^, S T ) T 
correspond to the three mutually perpendicular sides of a cube, and the neutrino basis 
L n = (S Ul , S U2 , S„ 3 ) T corresponds to another coordinate system (see Fig. [p. These two 
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FIG. 1: Geometric representation of the tribimaximal mixing in Eq. Q and 



coordinate systems are related to each other by rotations. One can perform two independent 
rotations to associate them. R\^{-k/A) and i?i2(sin~ 1 (l/y / 3)) transform the two independent 
bases into a common one. These practical operations are described below. R\ z {tx/A) mixes 
the second and third components of the basis L c and keeps the first one invariant to get a new 
basis (Si, S2, Ss) T , while i?i2(sin~ 1 (l/v / 3)) mixes the first and second components of L n , 
retaining the third one invariant to reach the same basis (Si, S2, S3) 71 . The mathematical 
expressions for relating (Si, S2, S3) 71 with the charged lepton basis (S e , S u , S T ) T and 
neutrino basis (S„ i; S„ 2 , S Ua ) T are shown as follows: 



S 2 



( 1 \ f ^\ 



4^/4) 



V S r y 



S 2 



i2i2(sm _1 (l/V3)) 



s, 



(8) 



Following the convention given in Ref. [3], when we discuss the geometry structure, we 
abbreviate the sides Si (S Vi ) as / (i/j) without causing any confusion. Here S 6|(UjT and S ul)V2tU3 
just refer to the corresponding geometrical quantities marked in Fig. [T]and are by no means 
the physical states. 

Comparing U tr iu = -R23( 7r /4)-Ri2(sin~ 1 (l/v / 3)) with Upmns = U\\J V , it appears that the 
two rotations R\ 3 (tt/4:) and i?i 2 (sin _1 (l/v / 3)) correspond to the mixing matrices for the 
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charged leptons and neutrinos, respectively. It is noted that in Eq. (jSJ), we only concern the 
mixing parts; thus, inserting 7°7 M (1 — 75) between ((Si, S2, S-s) T y and (Si, S2, S3) 71 which 
is irrelevant to our geometrical settings, we just derive the Lagrangian of weak interaction. 
We also would like to point out that in Eq. (jSJ), the high symmetry is assumed, and all 
quantities are indeed corresponding to the zeroth order ones [7], and then later when we 
introduce a deformation of the cube to break the tribimaximal symmetry, the concerned 
quantities would turn into the physical ones. 

Concretely, in Fig. (TJ sides OX, OY, and OZ represent e, /x, and r; and vi, u 2 , and 
correspond to OX' , OB, and OZ' , respectively. The \me-OZ' resides on the plane OZAY 
and spans an angle of 7r/4 with respect to the OZ axis, whereas line OX', line OZ', and 
line OB compose three-dimensional mutually perpendicular coordinate axes, and according 
the right-hand rule, we have an OB — OZ' — OX' system. The angle spanned between 
OX' and OX is 612- Therefore, the two rectangular coordinate systems transform from 
each other by two rotations about the axes OX and OZ', respectively. The right-handed 
rotation R\^(n/A) about the OX axis brings /i to OA and r to 1/3, and a second right-handed 
rotation i?i2(sin~ 1 (l/v^3)), with 612 = sin -1 (l/\/3), turns Vi into e and V2 into OA. Then, 
after performing the two successive operations, the basis (Si, S 2 , S3) 71 shown above can be 
directly read out as (e, OA, u 3 ) T ' . 

Although the tribimaximal mixing Ansatz is close to the experimental data and exhibits 
a striking symmetry, it is not the exact form of the PMNS matrix. Moreover, this symmetry 
demands #13 to be zero. If the tribimaximal symmetry is not exact, with the angles #23 and 
9 12 obviously deviate from the values determined by the symmetry, one has sufficient reason 
to believe that #13 should not be zero. In fact, the previous measurements set a lower bound 



for #13 as sin 2#i 3 < 0.15 [3[, and will be more precisely measured at the upcoming reactor 
experiments Daya Bay js| and Double Chooz {9]. 

It would be interesting to investigate how to break the tribimaximal symmetry from a 
theoretical aspect. Friedberg and Lee suggest to break the symmetry from the charged 



lepton side [10|, whereas He and his collaborators break the symmetry from the neutrino 
sector [ll I . Since the whole mixing matrix is a product of the two unitary matrices that, 
respectively, diagonalize the charged lepton and neutrino mass matrices as Upmns = U{U V , 
breaking from either side is just like climbing up Mount Everest from the south or north 
side as Lee comments 12[ . Their schemes to break the symmetry are algebraic. 
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FIG. 2: (color online) The three classes of symmetry axes for a cube, i.e., EE' , FF', and GG' . 

Instead, in this work, we propose to break the symmetry based on Friedberg and Lee's 
geometrical picture. Namely, we let the cube be slightly deformed and the nonzero 9 is value 
would emerge. Concretely, by deforming the geometric representation of the tribimaximal 
mixing, the angles would deviate from their ideal values; by fitting them to the data, we 
determine the deformation scale of the cube, and then by the new geometric shape #13 is no 
longer zero. 

The work is organized as follows. In Sec. II, we present our geometrical model of deform- 
ing the cube to get the 6*13 as a function of the other two mixing angles. Then, in Sec. Ill, 
we present our numerical results. The last section is devoted to our conclusion and some 
discussions. 



II. THE DEFORMED CUBE MODEL 

It is noticeable that the angle between lines OA and OB and that between OY and 
OA in the cube are precisely the two mixing angles of the tribimaximal matrix 6* 12 and 
6*23, respectively. For a perfect symmetry, which corresponds to a complete cube, we have 
9 12 = sin~ 1 (l/v / 3) and #23 = tt/4, which are determined by the geometry. It is then viable 
that a deformation would lead to the more realistic form of the PMNS matrix, and, thus, the 
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nonzero #13 would emerge. After this deformation, 812 and 823 are not the values given above 
anymore, but dependent on the form of the deformation. A cube is a kind of polyhedron 
with high symmetry described by a certain group, so a deformation of a cube should be 
regarded clS db symmetry breaking. 

Now, let us demonstrate how to deform the above cube. For choosing the deformation 
scheme, we set three principles: 

• There are three rotation axes for a cube as presented in Fig. [21 i.e., EE', EE', and 
GG'. Apparently, the axis GG' is related with the mixing angle 9 12 (/AOB), and 
the axis EE', which is parallel to the side OA, is related to #23 (/AOY). Then, the 
rest symmetry axis EE' may be related to the zero #13 in the tribimaximal mixing. 
Thus, after the supposed deformation, the three symmetries would be broken, and the 
value of the deformation angle is related to 9 13 . For simplicity, we just choose the 
deformation angle to be 6^3 . 

• For the tribimaximal mixing, #23 = vr/4 and #13 = 0, there exists the \i — r symmetry 
13] . The global fit [l^] gives #23 = 42.8°, which apparently breaks the \i — t symmetry. 



Thus, in the deformation, #23 (/AOY) should be changed from 7r/4 to some values in 
order to break the \i — r symmetry. 

• In Ref. [l^l, the global fits of # 12 and #23 are 34.4° and 42.8°, respectively. Thus, 
for the deformation, 9 12 (/AOB) and #23 (/.AOY) should be changed toward smaller 
values than sin _1 (l/v / 3) and 7r/4, respectively. 

Considering above principles, the simplest and most direct scheme to deform the cube 
is to slide the bottom face parallel to the top face, and a small angle would emerge, and 
this angle is identified as #13. The length of each side is unchanged during the slide. This 
operation is explicitly illustrated in Fig. |3] With the parallel slide, the bottom face becomes 
EEGH. To be consistent with Friedberg and Lee's picture and the principles we proposed 
above, we identify # 12 = /FAG, 9 2 3 = /BAF, and 9\% = /EAE . E is the point of 
intersection between side AE' and plane EFGH. AE\ and EqE\ are perpendicular to the 
diagonal line EG. 

Setting /E EG = a and in the rectangular triangle KtAAEiG, one has 

AG 2 = AE\ + E1G 2 = cos 2 13 + sin 2 13 sin 2 a + (V2 + sin 9 13 cos a) 2 . (9) 




FIG. 3: (color online) The sketch of the shift of the bottom face relative to the top face. 
In AAEF, 

AG 2 = 1+4 cos 2 #23 - 4 cos 6» 23 cos (vr - 9 12 - sin" 1 (2 sin 12 cos 6> 23 )) , (10) 
and in AEE F, 

E F 2 = E E 2 + EF 2 - 2E E ■ EF cos ZE EF (11) 

4cos 2 6> 23 - cos 2 #i 3 = sin 2 6>i 3 + 1 - 2sin# 13 cos(a + ^). (12) 

The geometrical relationship of the sides and angles in the deformed cube would determine 
Eq. Eq. QUI) , and Eq. (TP21 . From these equations, we can get an analytical expression 
of #i3, with respect to the other two mixing angles 9 12 and #23 as 

sin 2 #13 = 4 cos 4 #23 - 4 cos 2 #23 + 4 cos 2 #23 cos 2 (9 12 + sin" 1 (2 sin 9 12 cos #23)) + 1. (13) 

As we have mentioned above, a cube is a highly symmetric polyhedron that could be 
represented by a global 54 group This group has 24 elements classified in five conjugate 
classes. As shown in Fig. [21 a cube has three kinds of rotation axes, h = 2, h = 3, and 
h = 4, corresponding to FF', GG', and EE', respectively. All the rotation axes in the same 
h are equivalent. 
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FIG. 4: (color online) Comparison of our result (dashed line) with the Daya Bay expected sensitivity 
limit to sin 2 29±3 as a function of running time. 

It is notable that the angle between the axes of h = 2 and h — 4 is 7r/4, and the angle 
between the axes of h = 2 and h — 3 is sm~ 1 (yfT/3). In other words, the two angles are 
exactly that in the tribimaximal form of the PMNS matrix. Another angle corresponding 
to #i3 must be an angle between h = 4 and h — 4 itself, so #13 = 0. 

With the deformation, the symmetry of the cube is broken, and #13 acquires a nonzero 
value. We can then view i3 as the parameter representing the deviation from the cubic 
symmetry. It is then viable to define #13 as the angle between the "new" h = 4 axis and the 
"old" one. 

III. NUMERICAL RESULTS 

Ref. [3] presents the updated global fit to the three-generation neutrino mixing: 

012 = 34.4 ± 1.0°, 23 = 42.81^° ■ (14) 
Using the data as inputs, we obtain the numerical value of #13: 

sin 2 20i3 = 0.0238, 9 13 = 4.44°. (15) 
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FIG. 5: (color online) Comparison of our result (dashed line) with the Double Chooz expected 
sensitivity limit to sin 2 2#i3 as a function of running time. 

The errors of the fit would cause a theoretical uncertainty to #13: 

sin 2 2^3 = 0.0238jftK (16) 

The errors are rather large, and, in fact, to make sin 2 26^3 > 0, the lower bound shown in 
the above expression is set. This expression indicates that our prediction on 6* 13 is somehow 
sensitive to the input data and that in order to get more precise values of #13, more precise 
values of the input are needed. 

n n 

Two reactor neutrino experiments, Daya Bay |8[ and Double Chooz [9| , aiming to directly 
measure #13 are expected to reach a very high precision. We illustrate a relation of the 
expected sensitivity of the Daya Bay and Double Chooz as a function of the running time 
in Figs. H] and EE respectively, where we mark the central value of sin 2 2# 13 calculated in this 
work. 

Apparently, because of the high precision, the sin 2 29 ^ value from our model would be 
probed at the first run of the Daya Bay experiment. 

We show sin 2 2^13 as a function of #12 for #23 = 42.8° in Fig. El and the dependence on 
#23 for 6*12 = 34.4° in Fig. Particle Data Group (PDG) presents an upper bound of #13 as 
sin 2 2#i3 < 0.15 at CL=90% [sj, which is also marked in Figs. |6]and[71 From Figs. |6]and 
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FIG. 6: (color online) sin 2 2#i3 as a function of #12 for 623 = 42.8° from the toy model (solid line). 
The limit sin 2 26 13 < 0.15, CL=90% from PDG is plotted as the dashed line. 

[TJ we can see that the theoretically predicted value of sin 2 2# 13 is sensitive to both 12 and 
823. By the updated data, 812 and 823 are constrained within the range (31.9° — 36.5°) and 
(40.2° - 48.3°), respectively. 



IV. DISCUSSIONS AND CONCLUSIONS 



823 = 7r/4 and #13 = imply the so-called fi — r symmetry [13j embedding in the neutrino 
mass matrix, i.e., the mass matrix in the flavor basis has an obvious — v T permutation 
symmetry. This leads to the mass matrix with the form 



M 



( 



\ 



ABB 
BCD 
B D C 



\ 



(17) 



In Ref. 



11] , the authors discussed the soft breaking of the fi — T symmetry that arises from 



the Majorana mass term of the heavy right-handed neutrinos in the minimal seesaw model. 
From their fi — r symmetry breaking model, they derived a relation among the mixing angles 
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FIG. 7: (color online) sin 2 2#i3 as a function of #23 for 9\2 = 34.4° from the toy model (solid line). 
The limit sin 2 29 13 < 0.15, CL=90% from PDG is plotted as the dashed line. 



and Dirac CP phase 



7T 



023 r 



#13 COt #12 COS S. 



12 



(18) 



For the case that the Dirac CP phase 6 = and substituting the experimental fits #12 = 34.4° 
and #23 = 42.8° into Eq. (TTBl) . we obtain the value of #13 as 



sin 2 2#i 3 = 0.00276, #13 = 1.51°. 



(19) 



Instead, in a parallel work, Friedberg and Lee 10| suggested that one can break the \i — t 
symmetry at the charged lepton side in terms of a perturbation method, and they also 
showed that the breaking may lead to a nonzero #13. 

In this work, by deforming the cube that corresponds to a full tribimaximal form of the 
mixing matrix according the proposed principles, we derive the analytic relation among the 
three lepton mixing angles, and, taking the experimental data as inputs, we deduce the 
value of unknown mixing angle #13. The result gives sin 2 2#i3 = 0.0238, i.e., #13 = 4.44°. 
As noticed, our theoretical prediction favors smaller #13. As #13 is to be measured at the 
Double Chooz and Daya Bay experiments, our result indicates that in the future, there 
would be a great opportunity to fix the mysterious #13. The recent measurement of the 
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T2K collaboration 16[ indicates that sin 2 2#i3 falls in a rather wide range of 0.03(0.04) < 
sin 2 26*13 < 0.28(0.34), and the central value of our theoretical prediction is consistent with 
the lower bound set by the collaboration, and the error range is comparable. This value 
also does not contradict the new measurement by Main Injector Neutrino Oscillation Search 
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